Input Parameters Optimization in Swarm 
DS-CDMA Multiuser Detectors 

Taufik Abrao^, Leonardo D. Oliveira^, Bruno A. Angelico^, Paul Jean E. Jeszensky^ 

^Dept. of Electrical Engineering; State University of Londrina, Brazil; taufik@uel .br 
|http: / /www.uel . b r /pes soal/ taufik | www.uel.br/pessoal/taufik 
^Dept. of Telecomm. and Control Engineering, Escola Politecnica of the University of Sao Paulo, Brazil; 
' {leonardo, p j j}@lcs . poll . usp . br 

Q . ■^Federal University of Technology - Parana (UTFPR), Cornelio Procopio, Brazil; bangelicogutfpr . edu . br 

Q 



Abstract 

In this paper, the uplink dkect sequence code division multiple access (DS-CDMA) multiuser detection 
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Q , problem (MuD) is studied into heuristic perspective, named particle swarm optimization (PSO). Regarding different 

system improvements for future technologies, such as high-order modulation and diversity exploitation, a complete 
^ ■ parameter optimization procedure for the PSO applied to MuD problem is provided, which represents the major 

' contribution of this paper. Furthermore, the performance of the PSO-MuD is briefly analyzed via Monte-Carlo 



simulations. Simulation results show that, after convergence, the performance reached by the PSO-MuD is much 
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00 

. better than the conventional detector, and somewhat close to the single user bound (SuB). Rayleigh flat channel is 

initially considered, but the results are further extend to diversity (time and spatial) channels. 



^ ■ I. Introduction 

o3 I In a DS-CDMA system, a conventional detector by itself may not provide a desirable quality of service, 
once the system capacity is strongly affected by multiple access interference (MAI). The capacity of a 
DS-CDMA system in multipath channels is limited mainly by the MAI, self-interference (SI), near-far 
effect (NFR) and fading. The conventional receiver (Rake) explores the path diversity in order to reduce 
fading effect, but it is not able to mitigate neither the MAI nor the near-far effect [[T4ll . Il25l . In this context, 
multiuser detection emerged as a solution to overcome the MAI [|25l . The best performance is acquired 
by the optimum multiuser detection (OMuD), based on the log-likelihood function (LLF) ll25l . In ll24l 
it was demonstrated that multiuser detection problem results in a nondeterministic polynomial-time hard 
(NP-hard) problem. After the Verdu's revolutionary work, a great variety of suboptimal approaches have 
been proposed: from linear multiuser detectors [|25l . to heuristic multiuser detectors ||9l, fTl. 
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Alternatives to OMuD into the class of linear multiuser detectors include the Decorrelator [|23l . 
and MMSE [fT9l . Besides, the classic non-linear multiuser detectors include the subtractive interference 
cancellation (IC) MuD [flSl and zero-forcing decision feedback (ZFDF) [5]. In spite of the relatively low 
complexity, the drawback of (non-)linear, ZFDF, and hybrid cancelers sub-optimal MuDs is the failure in 
approaching the ML performance under realistic channel and system scenarios. More recently heuristic 
methods have been proposed for solving the MuD problem, obtaining near-ML performance at cost of 
polynomial computational complexity 0, [HI. Examples of heuristic multiuser detection (Heur-MuD) 
methods include: evolutionary programming (EP), specially the genetic algorithm (GA) ||4l, particle 
swarm optimization (PSO) [fT2ll . [|26ll . lfT6l and, sometimes included in this classification, the deterministic 
local search (LS) methods [[T3l . [[TtI . which has been shown to present an very attractive performance x 
complexity trade-off for low order modulations. 

Nevertheless, there are few works dealing with complex and realistic system configurations. High-order 
modulation Heur-MuD in SISO or MIMO systems were previously addressed in fllll, [l26l, [|l5l. In [|T5l . 
PSO was applied to near-optimum asynchronous DS-CDMA multiuser detection problem under 16— QAM 
modulation and SISO multipath channels. Previous results on literature [fT6l. [HI suggest that evolutionary 
algorithms and particle swarm optimization have similar performance, and that a simple local search 
heuristic optimization is enough to solve the MuD problem with low-order modulation [[TTI . However 
for high-order modulation formats, the LS-MuD does not achieve good performances due to a lack of 
search diversity, whereas the PSO-MuD has been shown to be more efficient for solving the optimization 
problem under M-QAM modulation [15]. 

Recent works applying PSO to MuD usually assumes conventional values for PSO input parameters, 
such [[TOl . or optimized values only for specific system and channel scenarios, such [fT6l for flat Rayleigh 
channel, [[TSl for multipath and high-order modulation, and [(H for multicarrier CDMA systems as well. In 
this paper, a wide analysis, with BPSK, QPSK and 16-QAM modulation schemes, and diversity exploration 
is carried out. 

This paper provides a quite complete parameter optimization of the PSO-MuD applied to DS-CDMA 
systems in Rayleigh channels with BPSK, QPSK and 16-QAM modulations. The text has the following 
organization: Section [Ilj presents the system model, including DS-CDMA, OMuD, and the PSO-MuD. 
The PSO parameter optimization is shown in Section [nil while Section [IVl exhibits some performance 
results in terms of Monte Carlo simulation (MCS). Finally, Section [Vl summarizes the main conclusions 
of this work. 
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II. System model 

In this Section, a single-cell asynchronous multiple access DS-CDMA system model is described for 
Rayleigh channels, considering different modulation schemes, such as binary/quadrature phase shift keying 
(BPSK/QPSK) and 16-quadrature amplitude modulation (16-QAM), and single or multiple antennas at 
the base station receiver. After describing the conventional detection approach with a maximum ratio 
combining (MRC) rule, the OMuD and the PSO-MuD are described. The model is generic enough to 
allow describing additive white Gaussian noise (AWGN) and Rayleigh flat channels, other modulation 
formats and single- antenna receiver. 



A. DS-CDMA 

The base-band transmitted signal of the fcth user is described as ||20l 

/ c °° 



Skit) 



^|§f:4^gk{t-^n 



(1) 



where Sk is the symbol energy, and T is the symbol duration. Each symbol d^^\ k = 1, . . . , K is taken 
independently and with equal probability from a complex alphabet set A of cardinality M = 2™ in a 
squared constellation, i.e., rf^*^ E A C C, where C is the set of complex numbers. Fig. [H shows the 
modulation formats considered, while Fig. |2] sketches the K base-band DS-CDMA transmitters. 
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Fig. 1. Three modulation formats with Gray mapping. 
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Fig. 2. Uphnk base-band DS-CDMA transmission model with K users. 
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The normalized spreading sequence for the k-th user is given by 

N-l 



= ^5^afcHp(t-nT,), 0<t<T, (2) 



n=0 

where ak{n) is a random sequence with N chips assuming the values {±1}, p{t) is the pulse shaping, 
assumed rectangular with unitary amplitude and duration Tc, with Tc being the chip interval. The processing 
gain is given hy N = T/Tc. 

The equivalent base-band received signal at gth receive antenna, g = l,2,...,Q, containing I symbols 
for each user in multipath fading channel can be expressed by 

I-l K L 
i=0 k=l 1=1 

/,»^^e^-^„M+^^(t), (3) 

with Ak = L being the number of channel paths, admitted equal for all K users, Tq^k,i is the total 

delajfl for the signal of the A;th user, £th path at gth receive antenna, e^'^? '"'^ is the respective received 
phase carrier; ?7^(t) is the additive white Gaussian noise with bilateral power spectral density equal to 
Nq/2, and h^q\(^ is the complex channel coefficient for the ith symbol, defined as 

_ ^V^) (A\ 
"'q,k,e — lq,k,e^ ' V'^-^ 



'■q,k,e ~ Jq,k,e^ 



where the gain 7^ [ ^ is a characterized by a Rayleigh distribution and the phase ^ ^ by the uniform 
distribution W[0, 27r]. 

Generally, a slow and frequency selective channel^ is assumed. The expression in ([3]) is quite general 
and includes some special and important cases: if Q = 1, a SISO system is obtained; if L = 1, the 
channel becomes non-selective (flat) Rayleigh; if h^^\ ^ = 1, it results in the AWGN channel; moreover, 
if Tg,fc,£ = 0, a synchronous DS-CDMA system is characterized. 

At the base station, the received signal is submitted to a matched filter bank (CD), with D < L branches 
(fingers) per antenna of each user. When D > 1, CD is known as Rake receiver. Assuming perfect phase 
estimation (carrier phase), after despreading the resultant signal is given by 



1 



(i+l)T 



ygl,i = ^ / rg{t)gk{t - Tg^k/)dt (5) 

JnT 

— ^kll'q^k,l^k '^^q,k,£~^ ^q,k,e~^ 'lq,k,£- 

The first term is the signal of interest, the second corresponds to the self-interference (SI), the third to 
the multiple- access interference (MAI) and the last one corresponds to the filtered AWGN. 



'Considering the asynchronism among the users and random delays for different paths 
^Slow channel: channel coefficients were admi 
{AB)c, the coherence bandwidth of the channel 



^Slow channel: channel coefficients were admitted constant along the symbol period T; and frequency selective condition is hold: 7^ >> 
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Considering a maximal ratio combining (MRC) rule with diversity order equal to DQ for each user, 
the M— level complex decision variable is given by 

Q D 

(6) 



q=i e=i 



k,i ■ Wq/fe,^, k = l,...,K 



where the MRC weights = 7gl £ e~■'^•'=■^ with ^^li and been a channel amplitude and phase 
estimation, respectively. 

After that, at each symbol interval, decisions are made on the in-phase and quadrature componentj^ 
of (^^ by scaling it into the constellation limits obtaining ^^'\ and choosing the complex symbol with 
minimum Euclidean distance regarding the scaled decision variable. Alternatively, this procedure can be 
replaced by separate v^M— level quantizers qtz acting on the in-phase and quadrature terms separately, 
such that 



rfl:)'^° = qtz(5R{e«n+jqtz(53U.. 



(7) 



•^rcal m 

for k = 1, . . . , K, and where ^reai and ^imag is the real and imaginary value sets, respectively, from the 
complex alphabet set A, and 3ft{-} and '^{■} representing the real and imaginary operators, respectively. 
Fig. [3] illustrates the general system structure. 
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Fig. 3. Uplink base-band DS-CDMA system model with Conventional receiver: K users transmitters, SIMO channel and conventional 
(Rake) receiver with Q multiple receive antennas. 



'Note that, for BPSK, only the in-phase term is presented. 
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B. Optimum Detection 

The OMuD estimates the symbols for all K users by choosing the symbol combination associated 
with the minimal distance metric among all possible symbol combinations in the M = 2™ constellation 
points IflSi . 

In the asynchronous multipath channel scenario considered in this paper, the one-shot asynchronous 
channel approach is adopted, where a configuration with K asynchronous users, / symbols and D branches 
is equivalent to a synchronous scenario with KID virtual users. 

Furthermore, in order to avoid handling complex- valued variables in high-order squared modulation 
formats, henceforward the alphabet set is re-arranged as ^rcai = ^imag = 3^ C Z of cardinality v^M, e.g., 
16-QAM (m = 4): d^^ ey = {±1, ±3}. 

The OMuD is based on the maximum likelihood criterion that chooses the vector of symbols d^, 
formally defined in (fT2l) . which maximizes the metric 



d°P* =arg max {l^fdjj, 



(8) 



where, in a SIMO channel, the single-objective function is generally written as a combination of the LLPs 
from all receive antennas, given by 

(9) 



9=1 

In the more general case considered here, i.e., K asynchronous users in a SIMO multipath Rayleigh 
channel with diversity D < L, the LLP can be defined as a decoupled optimization problem with only 
real-valued variables, such that 
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(11) 



KIDxl 



in Pq. (fTTI) is defined as 
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In addition, the G C^^^^^ is the despread signal in Eq. Q for a given q, in a vector notation, 
described as 



yyq,l,l yq,l,D) \yq,K,l yq,K,D) 



,(1) 



(1) 



,(1) 



{vqli ■ ■ ■ vIxd) ■■■(y, 



(I) 



q,K,l 



yq,K,D) 



(13) 



Matrices H and A are the coefficients and amplitudes diagonal matrices, and R represents the block- 
tridiagonal, block- Toeplitz cross -correlation matrix, composed by the sub-matrices R[l] and R[0], such 
that [l25l 

" R[0] R[l]^ ... 
R[l] R[0] R[1]T ... 
R[l] R[0] ... 



R 



... R[l] R[0] 
with R[0] and R[l] being KD matrices with elements 

1, if {k = u) and {i = I) 



(14) 



pIi,uP if < or 



</) 



P«,i,W if {k> u) OX {k = u4>l), 
0, 



if k>u 

1 

if k < u 

where a = (A; - !)£) + ^ 6 = (m - 1)D + / and fc, m = 1, 2, . . . , /sT 



a' 



(15) 



J = 1, 2 D; the cross-correlation 



element between the kth user, ith path and uth user, cith path, at gth receive antenna, pi ^ ^ ^, is 



Pk,e,u,d 



1 r 

^ Jo 



(16) 



The evaluation in ([8]) can be extended along the whole message, where all symbols of the transmitted 
vector for all K users are jointly detected (vector ML approach), or the decisions can be taken considering 
the optimal single symbol detection of all K multiuser signals (symbol ML approach). In the synchronous 
case, the symbol ML approach with / = 1 is considered, whereas in the asynchronous case the vector 
ML approach is adopted with / = 7 (/ must be, at least, equal to three (/ > 3)). 

The vector in (fTTI) belongs to a discrete set with size depending on M, K, I and D. Hence, the 
optimization problem posed by ([8]) can be solved directly using a m— dimensional (m = logg M) search 
method. Therefore, the associated combinatorial problem strictly requires an exhaustive search in A^^^ 
possibilities of d, or equivalently an exhaustive search in y'^^^^ possibilities of dp for the decoupled 
optimization problem with only real- valued variables. As a result, the maximum likelihood detector has a 
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complexity that increases exponentially with the modulation order, number of users, symbols and branches, 
becoming prohibitive even for moderate product values mK I D, i.e., even for a BPSK modulation format, 
medium system loading (K/N), small number of symbols / and D Rake fingers. 

C. Discrete Swarm Optimization Algorithm 

A discrete or, in several cases, binary PSO [fTTTl is considered in this paper. Such scheme is suitable 
to deal with digital information detection/decoding. Hence, binary PSO is adopted herein. The particle 
selection for evolving is based on the highest fitness values obtained through (flOl) and 

Accordingly, each candidate- vector defined like dj has its binary representation, bp[t], of size mKI, 
used for the velocity calculation, and the pth PSO particle position at instant (iteration) t is represented 
by the mKI x 1 binary vector 

b,[t] = [bjbj ■■■b; ■■■bf ]; (17) 

b; = ; e {o, i}, 

where each binary vector b^ is associated with one symbol in Eq. ([T2)) . Each particle has a velocity, 
which is calculated and updated according to 

vp[t + l]=u- vp[t] + 01 ■ Up,[t](bJ'=^'[t] - bp[t]) + 4>2 ■ Upjt](b5'=^'[t] - bp[t]), (18) 

where u is the inertial weight; Up^ [t] and Upj [t] are diagonal matrices with dimension mKI, whose 
elements are random variables with uniform distribution U E [0, 1]; b^'^'^'[t] and bp'^'*'[t] are the best global 
position and the best local positions found until the tth iteration, respectively; (pi and 02 are weight 
factors (acceleration coefficients) regarding the best individual and the best global positions influences in 
the velocity update, respectively. 

For MuD optimization with binary representation, each element in hp[t] in (fTSi) just assumes "0" or 
"1" values. Hence, a discrete mode for the position choice is carried out inserting a probabilistic decision 
step based on threshold, depending on the velocity. Several functions have this characteristic, such as the 
sigmoid function [11] 

s(v;jt]) = (19) 

where fp ;^[t] is the rth element of the pth particle velocity vector, = [f^^ ■ ■ -fp ^, ■ ■ -Vp^^], and the 
selection of the future particle position is obtained through the statement 

if u;^At]<s{v;^M, &;,.[t + i] = i; 

otherwise, 6^_^[t + 1] = 0, (20) 
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where 6p^[t] is an element of bp[t] (see Eq. ([TS])). and Up ^.[t] is a random variable with uniform distribution 

U G [0,1]. 

After obtaining a new particle position bp[t + 1], it is mapped back into its correspondent symbol vector 
dp[t + 1], and further in the real form dp[t + 1], for the evaluation of the objective function in 

In order to obtain further diversity for the search universe, the Knax factor is added to the PSO model, 
Eq. (fTSi) . being responsible for limiting the velocity in the range [±Vinax]- The insertion of this factor in 
the velocity calculation enables the algorithm to escape from possible local optima. The likelihood of a 
bit change increases as the particle velocity crosses the limits established by [±Vmax], as shown in Tab. 

m 

TABLE I 

Minimum bit change probability as a function of Vmax- 



^nax 


1 2 


3 


4 


5 


1 — 5'(V^nax) 


0.269 0.119 


0.047 


0.018 


0.007 



Population size V is typically in the range of 10 to 40 JH. However, based on lfT6l . it is set to 



P = 10 



0.3454 (^-^TiimKI - 1) + 2 



(21) 



Algorithm [T] describes the pseudo-code for the PSO implementation. 



III. PSO-MUD PARAMETERS OPTIMIZATION 

In this section, the PSO-MuD parameters optimization is carried out using Monte Carlo simulation. 
Such an optimization is directly related to the complexity x performance trade-off of the algorithm. A 
wide analysis with BPSK, QPSK and 16-QAM modulation schemes, and diversity exploration is carried 
out. 

A first analysis of the PSO parameters gives raise to the following behaviors: u is responsible for 
creating an inertia of the particles, inducing them to keep the movement towards the last directions of 
their velocities; (pi aims to guide the particles to each individual best position, inserting diversification 
in the search; 02 leads all particles towards the best global position, hence intensifying the search and 
reducing the convergence time; Vmax inserts perturbation limits in the movement of the particles, allowing 
more or less diversification in the algorithm. 

The optimization process for the initial velocity of the particles achieves similar results for three different 
conditions: null, random and CD output as initial velocity. Hence, it is adopted here, for simplicity, null 
initial velocity, i.e., v[0] = 0. 
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Algorithm 1 PSO Algorithm for the MuD Problem 



Input: d^D, -p, G, uj, 02, Kiax; Output: d^^° 
begin 

1. initialize first population: t = 0; 

B[0] = b*-"^ U B, where B contains {V — 1) particles randomly 
generated; 

h^''[0] = bp[0] and h^''[0] = b^O; 
Vp[0] = 0: null initial velocity; 

2. while t < G 

a. calculate ri(dp[t]), Vbp[t] e B[t] using 

b. update velocity Vp[t], p = 1, . . . through (fTsT i; 

c. update best positions: 
for p = 1, . . . , 7^ 

if l](dp[t]) > 17(4-'[t]), b^-'[t + 1] ^ bp[t] 
else b^'=^'[t + 1] ^ bj;"'[t] 
end 



if 3 bp[t] such that r2(d [t]) > f7(d^'=''[t]) 



A 



bb/»'[t + 1] ^ bp[t] 

else b^'=^'[t + 1] ^ ^''"[t] 

d. Evolve to a new swarm population B[t + 1], using ( l20b ; 

e. set t = t + 1. 
end 

3. bPSO =bb-'[G]; bPSO dPSO. 
end 



d^°: CD output. 

V: Population size. 

G: number of swarm iterations. 

For each [t] there is a bp [t] associated. 



In [fT6l, the best performance x complexity trade-off for BPSK PSO-MuD algorithm was obtained 
setting Vtnax = 4. Herein, simulations carried out varying V"max for different modulations and diversity 
exploration accomplish this value as a good alternative. This optimization process is quite similar for 
systems with QPSK and 16-QAM modulation formats. 

A. u Optimization 

It is worth noting that a relatively larger value for tu is helpful for global optimum, and lesser influenced 
by the best global and local positions, while a relatively smaller value for u is helpful for course 
convergence, i.e., smaller inertial weight encourages the local exploration [jH, [|2T1 as the particles are 
more attracted towards bp'^^'ft] and b^''''[t]. 

Fig. m shows the convergence of the PSO scheme for different values of to considering BPSK modulation 
and flat channel. It is evident that the best performance x complexity trade-off is accomplished with a; = 1. 



Many research papers have been proposed new strategies for PSO principle in order to improve its 
performance and reduce its complexity. For instance, in Q the authors have been discussed adaptive 
nonlinear inertia weight in order to improve PSO convergence. However, the current analysis indicates 
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that no further specialized strategy is necessary, since the conventional PSO works well to solve the MuD 
DS-CDMA problem in several practical scenarios. 

The optimization of the inertial weight, cu, achieves analogous results for QPSK and 16-QAM mod- 
ulation schemes, where to = 1 also achieves the best performance x complexity trade-off (results not 
shown). A special attention is given for (pi and 02 optimization in the next, since their values impact 
deeply in the PSO performance, also varying for each modulation. 



B. 01 and 02 Optimization 

1 ) BPSK Modulation: For Rayleigh channels, the performance improvement expected by 0i increment 
is not evident, and its value can be reduced without performance losses, as can be seen in Fig. |5l Therefore, 
a good choice seems to be 0i = 2, achieving a reasonable convergence rate. 
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Fig. 5. 01 optimization in Rayleigh flat channels with BPSK modulation, Eb/No — 22 dB, K = 15, (j)2 ~ 10, Knax = 4, and uj = 1. 



Fig. [61(a) illustrates different convergence performances achieved with 0i = 2 and 02 G [1; 15] for 
medium system loading and medium-high Eb/No. Even for high system loading, the PSO performance 
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is quite similar for different values of 02, as observed in Fig. [61(b). Hence, considering the performance 
X complexity trade-off, a reasonable choice for 02 under Rayleigh flat channels is 02 = 10. 
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Fig. 6. (/)2 optimization in Rayleigh flat channels with BPSK modulation, Vmax = 4, tj = 1, and 0i = 2; a) convergence performance with 
Ei/No = 22 dB and K = 15; b) average BER x K with Ei/No = 20 dB, G = 30 iterations. 

2j QPSK Modulation: Different results from BPSK are achieved when a QPSK modulation scheme 
is adopted. Note in Fig. |7] that low values of 02 and high values 0i delay the convergence, the inverse 
results in lack of diversity. Hence, the best performance x complexity is achieved with 0i = 02 = 4. 




Fig. 7. 01 and <j)2 optimization under flat Rayleigh channels for QPSK modulation, Eb/No = 22 dB, K = 15, uu = 1 and Vmax ~ 4. 



3) 16-QAM Modulation: Under 16-QAM modulation, the PSO-MuD requires more intensification, 
once the search becomes more complex due to each symbol maps to 4 bits. Fig. [8] shows the convergence 
curves for different values of 0i and 02, where it is clear that the performance gap is more evident with an 
increasing number of users and Eh/No. Analyzing this result, the chosen values are 0i = 6 and 02 = 1. 
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Fig. 8. 01 and <^2 optimization under flat Rayleigh channels for 16-QAM modulation, Et/No — 30 dB, K = 15, u = 1 and Vmax = 4. 

C. Diversity Exploration 

The best range for the acceleration coefficients under resolvable multipath channels (L > 2) for MuD 
SISO DS-CDMA problem seems (pi = 2 and 02 £ [12; 15], as indicated by the simulation results shown 
in Fig. |9l For medium system loading and SNR, Fig. |9] indicates that the best values for acceleration 
coefficients are 0i = 2 and 02 = 15, allowing the combination of fast convergence and near-optimum 
performance achievement. 

10-' 



10-' 

5 10 15 20 25 30 3 5 40 45 50 
Iterations 

Fig. 9. 01 and 02 optimization under Rayleigh channels with path diversity (L — D = 2) for BPSK modulation, Et/No = 22 dB, K — 15, 

OJ = 1, Kiax = 4. 




D. Optimized parameters for PSO-MuD 

As previously mentioned, the optimized input parameters for PSO-MuD vary regarding the system and 
channel scenario conditions. Monte-Carlo simulations exhibited in Section JV] adopt the values presented 
in Tab. |n] as the optimized input PSO parameters. Loading system C range indicates the boundaries for 
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TABLE II 

Optimized parameters for asynchronous PSO-MuD. 



Channel & Modulation 


£ range 


LO 01 02 Kiax 


Hat Rayleigh BPSK 


[0.16; 1.00] 


1 2 10 4 


Flat Rayleigh QPSK 


[0.16; 1.00] 


14 4 4 


Hat Rayleigh 16-QAM 


[0.03; 0.50] 


16 1 4 


Diversity Rayleigh BPSK 


[0.03; 0.50] 


1 2 15 4 



^ which the input PSO parameters optimization was carried out. For system operation characterized by 
spatial diversity {Q > 1 receive antennas), the PSO-lMuD behaviour, in terms of convergence speed and 
quality of solution, is very similar to that presented under multipath diversity. 

IV. NUMERICAL RESULTS WITH OPTIMIZED PARAMETERS 

In this section, numerical performance results are obtained using Monte-Carlo simulations. The results 
are compared with theoretical single-user bound (SuB), according to Appendix A, since the OMuD 
computational complexity results prohibitive. The adopted PSO-MuD parameters, as well as system and 
channel conditions employed in Monte Carlo simulations are summarized in Tab. Hill 

TABLE III 

System, channel and PSO-MuD parameters for fading channels performance analysis. 



Parameter 


Adopted Values 


DS-CDMA System 


# Rx antennas 


1, 2, 3 


Spreading Sequences 


Random, = 31 


modulation 


BPSK, QPSK and 16-QAM 


# mobile users 


K e [5; 31] 


Received SNR 


Eb/Nii e [0; 30] dB 


P50-MuD Parameters 


Population size, V 


Eq. (EB 


acceleration coefficients 


01 =2, 6; 02 = 1,10 


inertia weight 


UJ = 1 


Maximal velocity 


^nax 4 


Rayleigh Channel 


Channel state info. (CSI) 


perfectly known at Rx 




coefficient error estimates 


Number of paths 


L = 1, 2, 3 
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Fig. [To] presents the performance as a function of received E^/Nq for two different near-far ratio 
scenarios under flat Rayleigh channel. Fig. [ini(a) was obtained for perfect power control, whereas Fig. 
[TOl(b) was generated considering half users with NFR = +6 dB. Here, the BERAvg performance is 
calculated only for the weaker users. Note the performance of the PSO-MuD is almost constant despite 
of the NFR = +6 dB for half of the users, illustrating the robustness of the PSO-MuD against unbalanced 
powers in flat fading channels. 




Eb/No [dB] Et/No [dB] 

(a) (b) 
Fig. 10. Average BERAvg x Eb/No for flat Rayleigh channel with K = 15: (a) perfect power control; (b) NFR = +6 dB for 7 users. In 

scenario (b), the performance is calculated only for the weaker users. 



A. Diversity 

In the results presented here, two assumptions are considered when there are more than one antenna 
at receiver (Spatial Diversity): first, the average received power is equal for all antennas; and second, the 
SNR at the receiver input is defined as the received SNR per antenna. Therefore, there is a power gain 
of 3 dB when adopted Q = 2, 6 dB with Q = 3, and so on. The effect of increasing the number of 
receive antennas in the convergence curves is shown in Fig. [HI where PSO-MuD works on systems with 
Q = 1, 2 and 3 antennas. A delay in the PSO-MuD convergence is observed when more antennas are 
added to the receiver, caused by the larger gap that it has to surpass. Furthermore, PSO-MuD achieves 
the SuB performance for all the three cases. 

The exploitation of the path diversity also improves the system capacity. Fig. \TT\ shows the BERAvg 
convergence of PSO-MuD for different of paths, L = 1 2 and 3, when the detector explores fully the 
path diversity, i.e., the number of fingers of conventional detector is equal the number of copies of signal 
received, D = L. The power delay profile considered is exponential, with mean paths energy as shown 
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in Tab. |IV] [|20l . It is worth mentioning that the mean received energy is equal for the three conditions, 
i.e., the resultant improvement with increasing number of paths is due the diversity gain only. 
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Fig. II. Convergence performance of PSO-MuD, with K = 15, Eb/No = 15, BPSK modulation, (a) under asynchronous multipath slow 
Rayleigh channels and / = 3, for L = 1, 2 and 3 paths; and (b) synchronous flat Rayleigh channel, I — 1 and Q = 1, 2, 3 antennas. 



TABLE IV 

Three power-delay profiles for different Rayleigh fading channels used in Monte-Carlo simulations. 



Param. 


PD-l 


PD-2 


PD-3 


Path, i 


1 


1 2 


1 2 3 










T, 2T, 


nil] 


1.0000 


0.8320 0.1680 


0.8047 0.1625 0.0328 



Note there is a performance gain with the exploration of such diversity, verified in both the Rake 
receiver and PSO-MuD. The PSO-MuD performance is close to SuB in all cases, exhibiting its capability 
of exploring path diversity and dealing with SI as well. In addition, the convergence aspects are kept for 
all conditions. 

The PSO-MuD is also evaluated under channel error estimation, which are modeled through the 
continuous uniform distributions W [1 ± e] centralized on the true values of the coefficients, resulting 

7^ = U[l±e,]x 7« ; 9^^, = U[l±ee]x , (22) 

where and eg are the maximum module and phase normalized errors for the channel coefficients, 
respectively. For a low-moderate SNR and medium system loading {C = 15/31), Fig. [I2l shows the 
performance degradation of the PSO-MuD considering BPSK modulation, L = 1 and L = 2 paths or 
Q = 1 and Q = 2 antennas, with estimation errors of order of 10% or 25%, i.e., = eg = 0.10 
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Fig. 12. Performance of P^^-MuD with K = 15, BPSK modulation and error in the channel estimatidii,'' for (a) path diversity and (b) 
spatial diversity. 
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Fig. 13. Performance of Ps'tS-MuD with Et/No = 15 dB and BPSK modulation, for (a) path diversity ancl (b) spatial diversity. 



or = eo = 0.25, respectively. Note that PSO-MuD reaches the SuB in both conditions with perfect 
channel estimation, and the improvement is more evident when the diversity gain increases. However, 
note that, with spatial diversity, the gain is higher, since the average energy is equally distributed among 
antennas, while for path diversity is considered a realistic exponential power-delay profile. Although there 
is a general performance degradation when the error in channel coefficient estimation increases, PSO- 
MuD still achieves much better performance than the CD under any error estimation condition, being 
more evident for larger number of antennas. 

Fig. [13] shows the performance as function of number of users K. It is evident that the PSO-MuD 
performance is much superior then the CD scheme. 
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B. QPSK and 16-QAM Modulations 

Fig. [14] shows convergence comparison for three different modulations: (a) BPSK, (b) QPSK, and (c) 
16-QAM. It is worth mentioning, as presented in Tab. [Ill that the PSO-MuD optimized parameters is 
specific for each modulation. 
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Fig. 14. Convergence of PSO-MuD under flat Rayleigh channel, E^/No = 20 dB, and a) K 
K — 12 users with QPSK modulation and c) K — 6 users with 16-QAM modulation 



24 users with BPSK modulation, b) 



Similar results are obtained for E^/Nq curves in QPSK and 16-QAM cases. Nevertheless, Fig. [T5] shows 
that for 16-QAM modulation with (pi = 6, 02 = 1, the PSO-MuD performance degradation is quite slight 
in the range (0 < £ < 0.5), but the performance is hardly degraded in medium to high loading scenarios. 



V. CONCLUSION 

This paper provides an analysis of the PSO scheme applied to the multiuser DS-CDMA system, focusing 
on the parameters optimization of the algorithm. It was shown that to = 1 represents a good choice for 
the considered detection problem and configurations. 

Regarding the acceleration coefficients (0i and phi2) in Rayleigh flat channels, it was demonstrated 
that their choices depend on the modulation order. With BPSK (pi = 2 and (p2 = 10 represent a good 
choice. For QPSK 0i = 4 and 02 = 4 represented a good complexity x performance trade-off, while for 
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Fig. 15. PSO-MuD and CD performance degradation x system loading under 16-QAM modulation, in flat Rayleigh channel. 



16-QAM, it was observed that 0i = 6 and 02 = 1 provide a good result. However, in the latter case, the 
performance is not optimum for high system loading. With BPSK and Rayleigh diversity channels, it was 
shown that (pi = 2 and 02 = 12 to 15 provide a good convergence of the PSO. 

The PSO algorithm shows to be efficient for SISO/SIMO MuD asynchronous DS-CDMA problem 
when the input parameters are properly chosen. Under a variety of simulated/analyzed realistic scenarios, 
the performance achieved by PSO-MuD, except for high order modulation in the high system loading 
condition, was near-optimal. In the presence of channel errors, the PSO-MuD keeps much more efficient 
than conventional receiver with perfect channel estimation. In all evaluated system conditions, PSO-MuD 
resulted in small degradation performance if those errors are confined to 10% of the actual instantaneous 
values. 

Appendix 

A. Minimal Number of Trials and Single- User Performance 

The minimal number of trials {TR) evaluated in the each simulated point (SNR) was obtained based 
on the single-user bound (SuB) performance. Considering a confidence interval, and admitting that a non- 
spreading and a spreading systems have the same equivalent bandwidth {BW ^ ^ = -BPVspread ~ 
and thus, equivalently, both systems have the same channel response (delay spread, diversity order and 
so on), the SuB performance in both systems will be equivalent. So, the average symbol error rate for 
a single-user under A/-QAM DS-CDMA system and L Rayleigh fading path channels with exponential 
power-delay profile and maximum ratio combining reception is found in |l22l Eq. (9.26)] as 
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SERs,B = 2a -/?,)+ (23) 



a" 



1=1 

L 



y^PiPt X tan ^ ( ) - y^pi 

1^1 ^P'y . 



where: 

Pi= \ II 1 i - ^ I I , a 




a _ / ^i9qam 
Pi — \ ~^ — 5 9qam 



l + rn^Au ^■^'^^ 2(M-1)' 
and Vi = log2 M = mVl denotes the average received signal-noise ratio per symbol for the £th path, 
with V*^ being the correspondent SNR per bit per path. 

Once the lower bound is defined, the minimal number of trials can be defined as 



TR 



"errors 



SERsuB 

where the higher ?7^errors value, the more reliable will be the estimate of the SER obtained in MCS jHl. In 
this work, the minimum adopted ne^rors = 100, and considering a reliable interval of 95%, it is assured that 
the estimate SER C [0.823; 1.215] SER. Simulations were carried out using MATLAB v.7.3 plataform. 
The MathWorks, Inc. 
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